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Electron hopping models on the honeycomb lattice are studied. The lattice consists of two trian- 
gular sublattices, and it is non-Bravais. The dual space has non-trivial topology. The gauge fields 
of Bloch electrons have the U(l) symmetry and thus represent superconducting states in the dual 
space. Two quantized Abrikosov fluxes exist at the Dirac points and have fluxes 2n and — 2n, respec- 
tively. We define the non-Abelian 50(3) gauge theory in the extended 3d dual space and it is shown 
that a monopole and anti-monoplole solution is stable. The 5*0(3) gauge group is broken down to 
U(l) at the 2d boundary. The Abrikosov fluxes are related to quantized Hall conductance by the 
topological expression. Based on this, monopole confinement and deconfinement are discussed in 
relation to time reversal symmetry and QHE. 

The Jahn- Teller effect is briefly discussed. 

PACS numbers: 73.22.-f, 73.43.-f 
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I. INTRODUCTION 

There have been much scientific and technological in- 
terests in Carbon compounds. Notable examples are car- 
bon nanotubes and graphene, a single layer of carbon 
atomsi 2 -^. Integer quantum Hall effect has been reported 
in graphenei^. In both materials carbon atoms form the 
honeycomb lattice. 

The honeycomb lattice consists of two sublattices A 
and B which are the triangular lattice. A point in sub- 
lattice A and a point in sublattice B are not equivalent 
(see Fig. [JJ). Thus it is non-Bravais and the dual space 
has a quite non-trivial structure. 

In this paper we study the tight-binding model on the 
honeycomb lattice. The electronic states are Bloch waves 
which have a peculiar energy spectrum with the 2d Dirac 
zero modes and also the Id Dirac modes which have not 
been well-recognized. A U(l) gauge field is constructed 
from the Bloch waves. This kind of gauge field was 
originally introduced in 2c? periodic systems in a mag- 
netic field in the study of quantum Hall effcct(QHE) 4 . 
When time-reversal symmetry is broken, the gauge field 
is topologically non-trivial in the magnetic dual space 
and Chern numbers of the pertinent fiber bundle give 
quantized Hall conductance. 

In Sec. HH the tight-binding model is introduced which 
represents nearest-neighbor electron hopping on the hon- 
eycomb lattice. The Bloch states and the energy spec- 
trum are given. In Sec. IIII1 the dual space for the Bloch 
states are discussed. A convention to use two Brillouin 
zones is taken, and it is shown that it has a non-trivial 
topology being a torus T 2 . In Sec. IIVI the Dirac zero 
modes at E = are discussed. Especially it is shown 
that Bloch's theorem used in Sec. [Til is no longer valid 
and the states are high degenerate. These Dire modes 
are two-dimensional, but in Sec. [V] the existence of the 
Id Dirac modes is shown. The corresponding states are 
again highly degenerate and it is possible to construct 




FIG. 1: (Color online) The honeycomb lattice. Black and 
white circles are lattice points on sublattice A and B, re- 
spectively. Horizontal bonds between them are labeled by 
(n, m) and wave functions on the two ends are denopted by 
ipn, m (black circles) and (j) n ,m (white circles). The hopping in- 
tegrals for the horizontal bonds are t and those for the others 
are f . 



the strap dimer states which are strictly localized in one 
direction and extended in the other direction. 

In Sec. Ell the U(l) gauge field for Bloch electrons is 
introduced and then, in Sec. IVII1 superconductivity with 
quantized Abrikosov fluxes in the dual space is discussed. 
In Sec. IVIII| a massive case is considered in which there is 
no Dirac zero modes. Robustness of the fluxes is shown. 
In Sec. IIX1 the 2d dual space is extended to 3d In 
the U(l) gauge theory there is a magnetic flux line but 
there is no magnetic monopoles. On the other hand, 
there exists a monopole-antimonopole pair in the non- 
Abelian 50(3) theory which shows confinement. In this 
case, the 5*0(3) gauge group is broken down to U(l) at 
the 2d boundary. In Sec. |Xj the topological expression in 
QHE relates this monopole confinement and time reversal 
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symmetry. In Sec. IXI1 QHE of the honeycomb lattice 
is discussed. Energy spectra against magnetic fluxes per 
unit cell of the honeycomb lattice are shown for a number 
of anisotropic cases. They are compared with that of the 
square lattice. The extended 3d dual space is split into q 
subsapces and it is shown that monopoles are deconfincd 
due to breaking of time reversal symmetry. 

In the final Section lXIIl the Jahn- Teller effect is briefly 
discussed. 



II. TIGHT-BINDING MODEL ON THE 
HONEYCOMB LATTICE 
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The honeycomb lattice is an alternate lattice which 
consits of sublattice A and sublattice B. The two sublat- 
tices are triangular lattices (see Fig. [T]). Since a point in 
sublattice and a point in sublattice B are not equivalent 
the honeycomb lattice is not a Bravais lattice. Most of 
the nontrivial results which are discussed in this paper 
are immediate consequences of this fact. 

Denote wave functions on sublattices A and B as ip n , m 
and 4> n ,m, respectively, as shown in Fig. [TJ The tight- 
binding model with nearest neighbor hopping is 

— <fin+l,m-l — <f>n+l,m+l ~ ^n,m = Eljj n ^ m , 

— ^ra-ljm+l — 1pn-\,m-l — t^n,m — E(j) nmi (1) 

where hopping integrals of the horizontal bonds are t and 
those for the other bonds are 1. 

Sublattice A and sublattice B are the triangular Bra- 
vais lattices and each has a set of primitive vectors 



ai = (1,-1), 
a 2 = (1,1). 

Then the reciprocal lattice can be generated by 

b, = 7T(1,-1), 

b 2 = 7r(l,l), 



where 



bi ■ a,j — 2n5ij. 



(3) 



(4) 



Since the system is periodic, Bloch's theorem applies 
and the wave functions are written 



Vv, 



(5) 



Here the functions u(k) and v(k) are periodic in the dual 
space, namely 



u(k + R) = u(k), 
v(k +R) = v(k), 



(6) 



if R is a reciprocal lattice vector which is given by 

R = k 1 b 1 + k 2 b 2 , (7) 



FIG. 2: (Color online) The energy spectrum of the honeycomb 
lattice. 



with integers k\ and k 2 - 
Equations (fT]) and ((5} give 



" u(k) ' 


= E 


' u(k) ' 


v(k) 




v(k) 



h(k) 

where the 2x2 matrix hamiltonian is 
h(k) = 

with 

A(fc) 



A(fe) 
A*(fc) 



-t 



2e x cos k y . 



(2) From Eqs.© and © we have 



u(k) _ A(fc) 



and 



E = T \A(k)\. 



(8) 



(9) 



(10) 



(11) 



fl2) 



The energy spectrum for the isotropic case, t = 1, is 
shown in Fig. 2 and it is symmetric with respect to 
E = 0. This is the consequence of the fact that the 
tight-binding model |T]) is bipartite. 



III. TOPOLOGICAL STRUCTURE OF THE 
DUAL SPACE 

A metric of the lattice is not included in the last sec- 
tion. Thus k x and k y are dimensionless. Let us introduce 
a length scale a and scale k x and k y by, 

2 



A'. 



K„ 



V3a 
2 



- hi, 



(13) 



FIG. 3: (Color online) The first Brillouin zone of the honey- 
comb lattice and the Dirac zero modes. Dirac points K, K', 
and K" are equivalent since they differ by reciprocal lattice 
vectors each other. In the same manner Dirac points L, L', 
and L" are equivalent. 



In this scale the direct lattice has the hexagonal symme- 
try and the honeycomb lattice consists of two triangular 
sublattices whose lattice spacings are a. The primitive 
vectors of the reciprocal lattice ([3]) are scaled and given 
by 



W = ( 



2tt 2ti\ 



y/3a a 
2?r 2-k 



\/3a a 

The reciprocal lattice is generated as 
R! = kxW + k 2 b 2 ' 



(14) 



(15) 



where k\ and k 2 are integers. The first Brillouin zone 
is the Wigner-Seitz cell of the reciprocal lattice and it is 
shown in Fig. [3] 

A Brillouin zone is defined only for a Bravais lattice. 
So note that the first Brillouin zone here, in fact, is for 
sublattice A (or B) which is the triangular lattice. 

The corners of the first Brillouin zone are given by 



K': (0,- 



2tt 2vr 



47T ' 

3a * 



V : (0 



47T, 

3a ' 



K" 



2?r 



y/Za 3a 



■ ( 16 ) 



In Fig. [4] (a), two Brillouin zones are drawn. These two 
are equivalent, but this redundant representation is con- 
venient. This can be seen as follows: The vertical zone 
boundaries (green lines) of the first Brillouin zone are 
equivalent since they differ by a reciprocal lattice vec- 
tor (4-71-/^30,0). Therefore a horizontal line is a circle 
S 1 . On the other hand the zone boundaries of the first 
Brillouin zone connected by a vertical line are not equiv- 
alent. The vertical line, M, reaches an equivalent point 
only after going through the two Brilloin zones. 

The dual space is a torus, S 1 x M , which is shown in 
Fig. 2] (b). The top and the botom surfaces of S 1 x M 





(b) 



FIG. 4: (Color online) (a) The dual space which consists of 
two Brillouin zones. Green lines, red lines, and blue lines are 
equivalent, respectively for themselves. The line in ^-direction 
is topologically S 1 , and the line in y-direction M is also S 1 . 
(b) The dual space with non-trivial topology. 



have to be identified since they are related by a reciprocal 
lattice vector (0,47r/a). 

A direction in the dual space is equivalent to ones 
which are rotated by ±27r/3 due to the hexagonal sym- 
metry of the honeycomb lattice. Thus we call x-direction 
and the ones which are rotated by ±27r/3 as the S 1 di- 
rections. In a similar manner, y-direction and the ones 
which are rotated by ±27r/3 are called M directions. 



IV. DIRAC ZERO MODES 

At a glance a zero energy state seems to be given by 
A(fe) = from Eq. |jI5J). Note that, if so, h(k) in Eq. 
([9|) vanishes and the problem is ill-defined. 

Facing this situation let us go back to Eq. |T]) before 
applying Bloch's theorem §5§. For E — 0, tp n ,m and <\> n ^ m 
are decoupled and Bloch's theorem is no longer valid. 
Namely there is no phase relation between tp n ^ m and 4> n ,m 
because they are solutions of two independent equations, 

i>n-l,m+l + i>n-l,m-l + tl/>n,m = °, 

0. (17) 

Bloch's theorem is valid for each equation separately and 
Eq.([T7D is satisfied if 



0, cos k y 



or 



t 
t 



k x = ±7r, cos k y = - 

Note that there is no Dirac point if 
t > t* = 2. 



(18) 



(19) 



(20) 
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For the isotropic case, t = 1, we have 
k x = 0, ky = i-^-, 

7T 

k x = i7T, ky = zb — . 



(21) 



(22) 



These points correspond to the corner of the first Bril- 
louin zone (fT6)) (see Fig. [3]). 
Two degenerate states are 



and 



VVi.m (fe) 
<t>n,m{k) = 0, 

^n,m{k) = 0, 
fin,™ (fe) 



(23) 



(24) 



The first state is nonzero only on sublattice A and the 
second state is nonzero only on sublattice B. Since they 
are degenerate, any linear combination of the two states 
is an eigenstate, namely 



fn,m(k) = a(k)ip n<m (k) + P(k)<t> nm (k) 



(25) 



is an eingenstate for arbitarary functions a(k) and (3(k). 
Thus the states are highly degenerate. 

The key point is that non-Bloch forms may be taken 
for the eigenstates at the Dirac zero modes. 



V. ONE-DIMENSIONAL DIRAC MODES AND 
THE STRAP DIMER STATES 

If k y = ±7r/2, Eq. (flQ)) gives A = — t which is inde- 
pendent of k x . Instead of Eq.©, let us write 



Ipr, 



then Eq. (fTJ) reduces to 



= e ik » m i/} n 



"0 1" 




ipn 


= E 




1 




. ^™ _ 





Here n is arbitrary, and we have 

4>n = ±lpn 

E = Tt- 



(26) 



(27) 



(28) 



A state (4> n ,ipn) is shown in Fig. [5] and it is strictly 
localized in x-direction and extended in y-direction. 

States with different n's are degenerate, and therefore 
any linear combinations of them are also eigenstates. In 
particular, we choose ip n — e zk = cn ,(j) n — ±e lk * n in Eq. 
(1261) and 



,/, _ p ik x ii ik y m 

A , ik x n ik y m 

IJ n.m — - Lc ' c 




n-1 



n n+1 



FIG. 5: (Color online) A strap dimer state. It is only nonzero 
on the green circles. 



These staes are in a Bloch's form and extended in both 
x- and y-directions. From Eas- fTU]) and (fT2")) . Fermi ve- 
locities of these Id Dirac modes are given by 



Vf = ±2 cos k x 



(30) 



VI. GAUGE FIELDS FOR BLOCH ELECTRONS 

Gauge fields for Bloch electrons were constructed first 
in the context of quantum Hall effect where time reversal 
symmetry is broken by an external magnetic field^. Al- 
though time reversal symmetry is not broken here let us 
define, nonetheless, a field in the dual space by 



A(k) = -i[u(k)*V k u{k) + v{k)*V k v(k)}. 



(31) 



In order to clarify the gauge structure, consider overall 
phases of wave functions. If (u(k), v(k)) is an eigenvector 
of h(k), so is 



(u(ky lV {ky) = e if ^>(u(k),v{k)). 



= e i/(fe) 

The field (f3Tj) is transformed as 

A(ky = A(k)+W k f(k) 



(32) 
(33) 



Thus there is a formal equivalence between A(k) and the 
vector potential in electromagnetism. Therefore we call 
A(k) as a "gauge field" in the dual space. The gauge 
invariance comes from the freedom of choosing overall 
phases of wave functions. A "magnetic field" in the dual 
space is given by 



B{k) = V fc x A(k), 



(34) 



(29) 



which is gauge-invariant. 

At the Dirac points we have even stronger gauge in- 
variance. For a wave function in a form of (|25[) one may 
change phases of a(k) and (3(k) separately. Thus we 
have independent gauge symmetries on sublattice A and 
sublattice B. This phenomenon is caused by breaking of 
Bloch's theorem. 
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VII. SUPERCONDUCTIVITY AND 
QUANTIZED ABRIKOSOV FLUXES 



VIII. ABRIKOSOV FLUXES IN A MASSIVE 
CASE 



Take ratio of wave functions on sublattice A and sub- 
lattice B, 



z(k) 



(35) 



We have z(k) = u(k)/v(k) from the Bloch form of wave 
functions © and Eq. Ijlip gives 



z(k) = ± 



A(fc) 
|A(fe)f 



Thus write 



and 



u(k) = z(k)e 
w(Jfe) = e l3(k \ 



ig(k) 



Z{k) 



J8(k) 



(36) 



(37) 



(38) 



From Eq. (|37|). the gauge field l|3lj) is given by A(k) — 
—iz*(k)V 'kz(k) + 2Vkg(k). Since the second term is a 
pure gauge we can write, in a suitable gauge, 

A(k) = -iz*(k)V k z{k) = V fc 0(fc), (39) 

where Eq. (|38[) is used. This shows that gauge field 
A(k) is an exact 1-form in terns of 9(k) which are phase 
differences between wave functions on sublattiuce A and 
sublattice B. 

In this manner we have a U(l) gauge field theory and 
then it is natural to interpret it as a theory of supercon- 
ductivity. Let us identify A(fe) as the Ginzburg-Landau 
order parameter for superconductivity in the dual space. 
Then London current is j v (k) = —m 2 s A v (k), where m v 
is mass of the U(l) field. For the magnetic field in the 
dual space (f3"4"|) we have 

/ B(k)dS = [ V x A(k)dS 
Js Js 

A{k)di = I V fc 0(fc) = 0, (40) 
c Jc 

unless a coutour C contains a Dirac point. This is Meiss- 
ner effect] a magnetic field can not penetrate into bulk 
of a superconductor. 

If a contour C includes a Dirac point we have 



B{k)dS = 2vr, if C 3 K, 



/ B{k)dS = -2tt, if C 3 L, 
Js 



(41) 



for both the upper and the lower bands. Thus we have 
flux 2tt at K and flux — 2ir at iA These are quantized 
Abrikosov fluxes. The states are in Type II superconduct- 
ing with Abrikosov fluxes at the Dirac points. The total 
Abrikosov flux is zero for either the upper or the lower 
band. 



Let us include potentials m on sublattice A and —m 
on sublattice B. Hamiltonian ^ is modified as 



h(k) 



m A(Jfe) 
A*(fe) -m 



The eigenstates and energies are given by 
u(k) A(fe) 



v{k) 



z(k) 



E — m 



and 



E = ±^~r 



|A(fe)P 



(42) 



(43) 



(44) 



Since phases of z(k) remain the same, it can be shown 
that Eq. (|4Tj) is still valid. Thus we have quantized 
Abrikosov fluxes at the same points K and L with fluxes 
2tt and — 2-7T, respectively. 

Note, however, that there are no Dirac points at the 
Abrikosov fluxes. Therefore we conclude that the exis- 
tence of Abrikosov fluxes is nothing to do with Dirac zero 
modes. 



IX. MAGNETIC FLUXES AND MONOPOLES 

In Sec. IIII1 it is shown that the dual space is a torus 
T 2 . We extend it to a 3d space R 3 which is shown in Fig. 
[6] This 3c? space is periodic in z-direction with period 
47r I a. The 2d dual space is the surface of the cylinder 
with radius p = 2/^/ia. 

Let us define the Ginzburg-Landau Lagrangian in B?: 



L = J d 3 k{±F q 



aM 2 + -M 4 



-(45) 



where 



D lV a = (di-iAi) V a . 



(46) 



The Lagrangian (|45p is invariant under a U(l) gauge 
transformation: 



A, -> A, + d l A. 



(47) 



The potential term is chosen to give the field configura- 
tion in the 2d dual space which is described in Sec. IVIII 
In the 2d dual space there is a 2tt Abrikosov flux at K 
and a — 2-7T Abrikosov flux at L'. In the extended space, 
i? 3 , a flux line cannot have an end point due to topo- 
logical stability of the Abelian £7(1) gauge theory. Also 
a flux line costs an energy which is propotional to its 
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FIG. 6: (Color online) The extended 3d dual space R 3 — 
(p,9,z). A 2ir Abrikosov flux is at K and a —2tt Abrikosov 
flux is at L. Points K' and L' are equivalent to K and L, 
respectively. A monopole is at the blue circle and an anti- 
monopole is at the red circle in the non-Abelian SO(3) gauge 
theory. The 3d dual space R 3 is periodic in z-direction with 
period 4tv /a. 

length. From these it can be seen that a magnetic flux 
line is a straight line which goes through from L' to K. 
(There is another flux line which goes through from L to 
K', but this is equivalent to the former.) Thus there is 
no magnetic monopole. 

A much more interesting subject is the non-Abelian 
gauge theory. Let us consider the 50 (3) gauge symmetry 
which has U(l) as a subgroup. 

L = J d 3 k{±(F tj y + l{D iV f + V(<p)], (48) 

where 

Ffj = d l A a J ^d J A b l -e abc A\A c v 

f = (v 1 , f 2 , <P 3 ), 
Ditp a = d t tp a -e abc A b <p c . (49) 

The potential V(ip) is chosen to give the U(l) field con- 
figuration on the 2d dual space which is described in Sec. 
IVIII Thus the 5*0(3) gauge symmetry is broken down to 
the subgroup U(l) on the surface. 

This symmetry breaking is forced on the surface and 
is not spontaneous as in the spontaneously broken gauge 
theory of 't Hooft^ and Polyakov. The U(l) gauge struc- 
ture of Bloch electrons on the honeycomb lattice deter- 
mine the symmetry on the 2c? surface and they break the 
50(3) gauge symmetry. 

As shown by 't Hooft 9 monopoles are allowed in the 
50(3) gauge theory. From the symmetry of the dual 
space R 3 we conclude that there exists a pair of monopole 
and antimonopole as shown in Fig. [5] 



X. HALL CONDUCTANCE AND MONOPOLE 
CONFINEMENT 

The Hall conductance in units of e 2 /h from a single 
band is given by the topological expression: 

<? xy = ^ J d 2 k[V k x A(k)] z , (50) 

where the integral is over the first Brillouin zone-K This 
is always an integer and called a TKNN integer—. In the 
present problem we do not expect that the Hall effect 
takes place because there is no magnetic field and time 
reversal symmetry is not broken. 

On the other hand, as shown above, there exists a pair 
of monople and antimonoplc. There is neither an iso- 
lated monopole nor an isolated antimonopole. This fact 
and the vanishing of the Hall conductance is related by 
the topological expression of Hall conductance (f50|) . The 
R.H.S represents the total flux (sum of the fluxes at Dirac 
points K and L) through the 2d dual space. 

Therefore we conclude that the monopole and anti- 
monopole pair confinement is forced by time reversal sym- 
metry. 

XI. QUANTUM HALL EFFECT AND 
MONOPOLE DECONFINEMENT 

In order to describe electron hopping on the honey- 
comb lattice in a magnetic field, phase factors are intro- 
duced in Eq. fTJ , 

— 4>n+l,m-l — 4>n+l,m+l — te l7r$m ' (j) n ,m = Elp n , m , 
—lpn-l,m+l - lpn-l,m-l ~ ie i,r$m V>n,m = E(j> n>m .(51) 

Here a magnetic flux through a unit hexagon is given by 
27r<I>. The energy spectra are shown in Fig. [Jj for (a) 
t = 1.0, (b) t = 1.4, (c) t = 2.0, and (d) t = 2.4. If 
there is no magnetic field, $ = 0, the spectra consist of 
an upper band and a lower band. As shown in Sec. IIV( 
the two bands touch at the Dirac points for t < t c = 2. 
This situation is similar to the case of the square lattice 
with $ = 1/2 (see Fig. EJ). For t > t c — 2 there is no 
Dirac point but there is a gap around E = (see Fig. [7] 
(d))- 

If <!> = p/q, where p and q are mutually prime integers, 
the translational symmetry of the honeycomb lattice in 
j/-direction is broken. Equation (|5ip becomes periodic 
with period q in y-direction. As a consequence the 2d 
dual space is spilt into q subspaces and both the upper 
and the lower bands are split into q subbands (see Fig. [7]). 
This phenomena has been known for the square lattice 
(Fig. [H]) and let us call it as the Hofstadter mechanism 
of band splittings^. 

When the Fermi energy is in a gap, the Hall conduc- 
tance is given by a sum of contributions from the sub- 
bands below the Fermi energy. A contribution from a 
single subband is given by the topological expression (|50|) 
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FIG. 7: (Color online) The energy spectra of the honeycomb 
lattice in magnetic fields for (a) t = 1.0, (b) t = 1.4, (c) 
t — 2.0, and (d) t = 2.4. A magnetic flux through a unit cell 
is given by 27r$. 



and thus quantized. Quantized Hall effect of the honey- 
comb lattice following these topological considerations is 
discussed in RefA 

If t > t c = 2, there is a gap around E — with zero 
Hall conductance (see Fig. [Tfd)) . In fact, both the up- 
per and the lower parts of the spectrum consist of q sub- 




0.2 0.4 0.6 0.8 1 
$ 

FIG. 8: (Color online) The energy spectrum of the square 
lattice in a magnetic field. 



n'+l,m'+l 




FIG. 9: (Color online) The square lattice for bonding states 
and ant-bonding states. 



bands and they are topologically equivalent, respectively, 
to that of the square lattice which is shown in Fig. [5] In 
order to see this, first consider the case <£> = 0. Take sum 
and difference of the two equations in Eq. ([I]) , 

4>n+l,m-l + <Pn+l,m+l + Wn—X,m+1 + V'n-l.m-l 
= ~{E + t)(t/} n>m +</>„, m ), 

4>n+l,m-l + 4>n+l,m+l ~ ^n-l.-m+l — V'n-l.m-l 

= -(E-t)(iP n>m -<j> ntm ). (52) 

For J > 1 we have bonding states, 4>n,m — ip n .m, and 
anti-bonding states, <j) n _ m — —ip n .m- For those states 
Eq. ([52]) represents a tight-binding model on the square 
lattice which is shown in Fig. [51 The solutions for those 
two sectors are: for bonding states, 



lb , i = 6 i / — Jkm'+ihim! 
Yn ,m' — Yn ,m' — c ; 

E = — t — (cos k\ -\- cos h>2 ) ) 
and for anti-bonding states, 

,/, , , — _J> , , — ikxn +ik 2 m 
Yrv ,m' — Yn ,m' — c j 

E = t — (cosfci + COSfo), 



(53) 



(54) 



where (n' , to') are the coordinates of the square lattice in 
Fig. [5]which are rotated from those in Fig. [TJ and k\ and 



ki are the corresponding k vectors which are in region: 

— 7T < ki < 7T, 
-7T < fc 2 < 7T. (55) 

Let us apply a magnetic field which gives flux 2ir& = 
2irp/q per unit cell, where p and q are mutually prime 
integers. Then each dual space is spilt into q subspaces; 
we have q subbands for each sector. In Fig. [7] (d) the 
lower Hofstadter butterfly is in the bonding sector and 
the upper Hofstadter butterfly is in the anti-bonding sec- 
tor. 

The quantized values of the Hall conductance for each 
sector are given as follows^: If a Fermi energy is in 
rth gap from the bottom, integer Hall conductance t r is 
determined by the Diophantine equation 

pt r + qs r = r, (56) 

where the integers must satisfy 

< r < q. (57) 

The integer Hall conductance carried by rth band from 
the bottom is given by a xy = 2n(t r — t r _i) and it is re- 
lated to the total flux of the dual subspace by the topolog- 
ical relation ([511)1 . Thus Hall conductances of subbands 
in a magnetic field are nonzero as it should be, because 
time reversal symmetry is broken. 

In SeclIX| the 2d dual space of the honeycomb lattice is 
extended to 3d space R 3 . The existence of a monoplole- 
antimonopole pair in R 3 is shown in the SO (3) gauge 
theory. In a magnetic field, the 2c? dual space is split into 
q subspaces by the Hofstadter mechanism. Accordingly 
R 3 is also split into q subspaces R 3 . The total Abrikosov 
flux of a single subband is nonzero since the correspond- 
ing quantized Hall conductance is nonzero. This implies 
that numbers of monopoles and anti-monopoles are un- 
equal in R 3 . This is monopole-antimonopole deconfine- 
ment due to time reversal symmetry breaking. 

Note that sum of Hall conductances of all the subbands 
vanishes. This can be seen from the fact that the only 
solution of the Diophantainc equation (|56|) with Eq. (|57|) 
is t r = for r = q. This implies that the total number of 
monopoles and antimonopoles are equal in q subspaces 
R 3, s. Thus, when a magnetic field is applied, monopoles 
and antimonopoles are created in pairs. But they are 
unpaired and distributed to subspaces R 3, s. 

XII. THE JAHN-TELLER EFFECT 

The Jahn- Teller effect was first predicted as a very gen- 
eral phenomenon 11 . It is the intrinsic instability of an 
electronically degenerate compound against lattice dis- 
tortions that remove degeneracy. Thus a dynamical lat- 
tice symmetry breaking takes place to a lower lattice sym- 
metry. 
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FIG. 10: The honeycomb lattice with a period 3 modulation 
of hopping integrals with tA, is, and tc in y-direction. The 
hopping integrals for the horizontal bonds are 1. 



As discussed in Sec. HVl there arc two Dirac modes 
with E = in the electron hopping model on the honey- 
comb lattice. Thus if the Fermi energy is at E = 0, the 
ground state is fourfold degenerate. (Each Dirac mode 
is doubly degenerate.) The condition for the Jan- Teller 
effect is satisfied but the whole problem is a rather dif- 
ficult one since all the interactions involoving electrons 
and nucleus have to be taken into account properly. 

Instead of doing it, let us focus on a simpler problem 
of considering ground state energes of electronic systems 
with certain periodic modulations of hopping integrals. 
From Eq. (fT6|) the vector connecting the Fermi surfaces 
(Dirac points K and L) is given by 

47T 

kper3 = (0, — ). (58) 

Then we would expect Fermi surface instability modes 
in M directions with k per ^ whose magnitude is a third of 
the period of the dual space(see Fig. 0]). 

Gap openings for period 3 modulations in M directions 
are numerically examined (see Fig. I10p . Examples are 
shown in Fig. 11 and we have numerical evidences of gap 
openings with period 3 modulations. 

The energy spectra for period 2 modulations of hop- 
ping integrals can be obtained analytically and it is con- 
firmed that there is no gap opening for modulations in 
neither the S 1 directions nor M directions. Periods 4 and 
5 modulations are examined numerically and there are 
no gap openings. From these investigations we conclude 
that gap openings take place only period 3 modulations 
of hopping integrals in M directions. 

The situation described above has a certain similarity 
to the Peirerls instability of a one-dimensional half-filled 
band in which 2kp instability leads to gap openings and a 
dimerization of the lattice^. Note that density of states 
for the Id lattice is finite at Ep = 0, but it vanishes for 
the honeycomb lattice. In this respect the honeycomb 
lattice is easier to be distorted compared with the 1-d 
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FIG. 11: Examples of gap openings with period 3 modulations 
of hopping integrals, (a) tA = tc — 1-02 and ts = 0.98. (b) 
t A = to = 1-05 and t B = 0.95. 



lattices. 

In Sec. [V]the existence of Id Dirac modes is explicitly 
demonstrated. Therefore from these modes, we expect 
to have another types of Jahn- Teller effect if Ep = ±t- 

The Id Dirac modes are on lines k y = ±7r/2, and if 
Ep = ±i, the two Fermi surfaces nest with a nesting 
vector 



2tt 

finest — ("j J) 

a 



(59) 



(it is scaled by Eq. (|13[) L Thus a possibility of charge 
density wave(CDW) or spin density wave(SDW) with pe- 
riod 2 in M direction has to be examined. 



Acknowledgments 

We are especially grateful to J. Friedel for pointing 
out the existence of Jahn- Teller effect in graphen, and 
to T. Eguchi for informing the author about the work of 
't Hooft on magnetic monopoles. We thank D. Tobe for 
help with computer calculations. 



1 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. 
I. Katsnelson, I. V. Grigorieva, S. V. Dubonos and A. A. 
Firsov, Nature 438, 197 (2005) 

2 Y. Zhang, Y.-W Tan, H. Stormer and P. Kim, Nature 438, 
201 (2005). 

3 C. Berger, Z. Song, X. Li, X. Wu, N. Brown, C. Naud, D. 
Mayou, T. Li, J. Hass, A.N. Marchenko, E. H. Conrad, P. 
N. First, W. A. de Heer, Science, 312, 1191 (2006). 

4 M. Kohmoto, Ann. of Phys. 160, 343 (1985). 

5 In physical units, 2n flux corresponds to the flux quantum 
$o = hc/e (or <I>o = hc/2e for superconductor). 



6 M. Kohmoto, Phys. Rev. B39, 11943 (1989). 

7 D. Thouless, M. Kohmoto, P. Nightingale, and M. den Nijs, 
Phys. Rev. Lett. 49, 405 (1982). 

8 Y. Hasegawa and M. Kohmoto, Phys. Rev. B74, 155415 
(2006). 

9 G. 't Hooft, Nucl. Phys. B79, 276 (1974). 

10 D. Hofstadter, Phys. Rev. B14, 2239 (1976). 

11 H.A. Jahn and E. Teller, Proc. Roy. Soc. 161, 220 (1937). 

12 R.E. Peierls, Quantum Theory of Solids (Oxford: Claren- 
don Press, 1956). 



